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ABSTRACT 

Even if the luminosity distance as a function of redshift is obtained accurately using, 
for example, Type la supernovae, the equation of state of the Universe cannot be 
determined uniquely but depends on one free parameter Q^o = /c/(ag-ffg), where ao 
and Hq are the present scale factor and the Hubble parameter, respectively. This 
degeneracy might be resolved if, for example, the time variations of the redshift of 
quasars are measured as proposed recently by Loeb. Therefore the equation of state 
of the Universe (or the metric of the universe) might be determined without any 
theoretical assumption on the matter content of the Universe in future. 
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1 INTRODUCTION 

To determine the structure and dynamics of an astrophysical system, the equation of state is usually necessary. For example, 
consider the structure of a spherical neutron star. If the pressure p is known as a function of the density p, we can determine the 
gravitational mass M and the radius R of the star as a function of the central density p c by solving the Oppenheimer-Volkoff 
equation (Oppenheimer & Volkoff 193!: ) . This means that we can determine the mass-radius relation M(R) theoretically in 
principle. However, the equation of state relevant to the neutron star is not established yet, although it may be determined by 
Quantum Chromo Dynamics in future. Therefore the mass-radius relation of the neutron star is not known well theoretically 
at present. 

Observationally, however, the mass-radius relation of neutron stars may be determined in the near future. If gravitational 
waves from the coalescing binary neutron stars are detected by the LIGO/VIRGO/GEO/TAMA network which will be in 
operation around 2000, the mass of each neutron star as well as its radius may be determined by analyzing the waveform in 



the last three minutes of the binary inspiral (rhornel995). In general the mass of each observed neutron star can be different 
so that there is a chance to determine M(R) observationally. In this case, as shown by Lindblom (1991), the equation of state 



of the high density matter can be determined from M (R) . This is, in a sense, the inverse problem. 

Now let us consider an isotropic and homogeneous universe which describes the Universe quite well in a global sense. 
The amount of the radiation p r (z) in the Universe is well known as a function of the redshift z from the present temperature 
of the cosmic background radiation. The amount of baryons pb(z) has a constraint from the big bang nucleosynthesis. We 
know that dark matter should exist but we do not know well its density Pd(z) or pressure Pd(z)- A cosmological constant 
may exist but introducing a non-zero cosmological constant needs a fine-tuning of the vacuum energy, and at present we do 
not have any convincing explanation for why such an extremely small value of the cosmological constant (in Planck units) 
is required. Since we do not know the equation of state of the matter in the standard model well, we cannot determine the 
scale factor as a function of time theoretically. Moreover recently several authors have considered the more general equation 



of state for a dark component called x-matter or quintessence and have explored its cosmological implications (Turner & 



Whitcl997 ; Chiba, Sugiyama fc Nakamural997 Chiba, Sugiyama & Nakamural998; Caldwell, Dave & Steinhardtl998; Chiba 
jr^akamural99c| ) . The situation is worse than the neutron star case; the equation of state in the expanding universe is almost 
unknown theoretically. 

Observationally we have several quantities such as the luminosity distance, the angular diameter distance and number 
counts as a function of t he redshift. Among these the luminosity d i stance cIl(z) m a y be determined q uite accurately by 
using Type la supernovae ( Garnavich et al.1998 Perlmutter et al.l99cf Riess et al.199 |; ISchmidt et al,1998| ). In the future the 
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afterglows of a certain class of gamma ray bursts might serve as a standard candle( Cohen & Piranl997). Therefore in this 
paper we assume that a quite accurate luminosity distance may be obtained and examine whether the equation of state of the 
expanding universe can be determined uniquely. Namely we discuss the inverse problem in cosmology. By "equation of state 
of the universe", we mean the relation between the total energy density of cosmic matter and the total pressure. Nearly three 
decades ago, Weinberg studied the possibility of determining the metric from the observed luminosity distance, with negative 
conclusion ( |Weinbergl970| ). We shall argue that the improvement of observational techniques now enables us to determine 
the metric of the universe (or the equation of the state of the universe) directly from the observational data. 

In section 2 we show that the equation of state can be determined if the scale factor is given as a function of time. We 
also show that the scalar field potential can be determined similarly. In section 3 we discuss how to determine the metric of 
the universe from the luminosity distance as a function of the redshift z. 



2 MATTER FIELD IN TERMS OF SCALE FACTOR 

Consider any given scale factor a(t) which is a monotonically increasing function of time. We denote the inverse function of 
a(t) as t(a). Then every function of t can be considered as a function of a. For example the Hubble parameter H can be 
written as 

:-(-=r-*w (1 > 

The metric of our isotropic and homogeneous universe is given by 

ds 2 = -dt 2 + a{t) 2 {d X 2 + f(xfdn 2 ) (2) 

where f(x) = X, sm h(x) an d sin(x) for flat, open and closed universes, respectively. From the general classification of energy 
momentum tensors (Hawking & Ellis 1973), the form of the energy momentum tensor compatible with the metric (2) is of 
Type I with pressure p\ = p2 — pz — p(a). We express the total energy density as p(a). Then p can be expressed by p as p(p) 
which we call the equation of state of the Universe. The Einstein equations are 



H(af = •fp-^, (3) 

H(a)+H(af = -£( p + 3p), (4) 

p = -3H(p+p), (5) 

where k 2 — 8ttG and k — 0, —1 and 1 for flat, open and closed universes, respectively. First, from Eq.(^) p as a function of a 
can be written as 

K 2 p(a) = 3 (H{a) 2 + ±) . (6) 
Then, from Eq.(^) p as a function of a can also be written as 

K 2 p(a) = - faff (a) + 3H(a) 2 + . (7) 



k k 



Since p(a) and p(a) are given, we can determine the equation of state p(p) . Note that we can consider even the case where 
the weak energy condition (Hawking & Ellis 1973) is violated: p + p<0orp<0. 

Next consider the minimally coupled scalar field (f> with the potential V(4>), which is an example of the general energy 
momentum tensor described above. Then 

P =±4, a + v(4>), (8) 

p = \tf-V{4>), (9) 

4> + 3>H^ = -V'. (10) 

In this case the weak energy condition is satisfied. We assume 4> > for simplicity though it is not necessary. From Eqs.Q 
and (^) we immediately have, 

k 2 V=(3H 2 +H+^, (11) 

«v-'(2".)*~ a (*-£)■ p» 

From Eq.@ 

V can be written as a function of a. Also, from Eq.@ cf> can be written as a function of a. Therefore, V can 
be written as a function of d>. 
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DETERMINING THE METRIC FROM OBSERVATIONS 



3.1 Matter field in terms of the luminosity distance 

From observations we might know, for example, the luminosity distance as a function of redshift cIl(z) accurately by using 



Type la supernovae( Garnavich et al.1998 



Perlmutter et al.1998 



Riess et al.1998 



Schmidt et al.1998). In this section we 



therefore will regard the Hubble parameter H(z) as a function of the redshift instead of as a function of the scale factor in 
the previous section. p(z) and p(z) are then given by 



K 2 p(z) = 3 (H(z) 2 - 
K 2 p(z) = -3H(z) 2 



(1 

2(1 + z)H(z 



z) H £l k0 ) , 
dH 

dz 



(l + z) 2 Hfa k0 , 



n 2 H 2 ' 
a a n 



(13) 
(14) 

(15) 



where ao and Ho are the present scale factor and the Hubble parameter, respectively. It is apparent that in general the 
equation of state of the Universe depends on Qko. To determine H(z) we will use the luminosity distance. Since all cosmological 
observations are made on the past light cone, the argument is similar for other distance indicators. The luminosity distance 
rfi(«) is defined by 



d L (z) 



1 

X = — 



ao(l + *)/(*) 
da 



(1 



/„ H(u) 

Then H(z) can be written in terms of r(z) as 



H(z) = (dr/dz)- 1 y/l - r(z) 2 H 2 Q. 



HI- 



(16) 
(17) 

(18) 



Note here that the above formula is valid irrespective of the sign of k. Since r(0) = 0, Ho can be determined irrespective 
of fifeo- One may think that fifco can be determined only from r(z). However, this is not the case. To show this, we rewrite 
Eq.dl8l) as 



dr 



r(z) 2 H 2 Q k0 . 

From the second derivative of Eq.(|l£ 
d 3 r 



(19) 



r(0) 



(20) 



we have 

dz 3y "' dz 2 \HJo Ho ' 

Eq.(|2o| shows that we cannot determine Qko without the knowledge of the second derivative of the Hubble parameter. This is 
due to the fact that as far as the expansion of the universe is concerned, the effect of the curvature is equivalent to "matter" 
with the equation of state p — —p/S. 

Now using r(z) we express p{z) and p(z) explicitly as 



k p(z) — 3 

K 2 p(z) = — 



(dr/dz) 2 
3 



+ (1 + *) 



* + < 1 + *>5 



(dr/dz) 
1 

(dr/dz) 2 



HoD.ko 



(1+*) 



3r 2 



/, \ d 
dz 



(dr/dz 



(21) 



(22) 



(dr/dz) 2 ' v ' ' dz \(dr/dz) 2 J ' ' (dr/dz) 2 

The above equations show that the equation of state of the Universe cannot be determined uniquely from the luminosity 
distance di,(z) = (l + z)r(z) but depends on one free parameter flko- This degeneracy was first pointed out by Weinberg ( 1970). 



This means that even if the luminosity distance is observed accurately as a function of z, in order to determine the equation of 
the state of the Universe (or the metric) some assumption about Q,ko is needed. Intuitively, this can be understood as follows. 
Since the luminosity distance or other cosmological distance measure only gives information on the past light cone, it is not 
enough to infer what happens inside the light cone. 



3.2 Breaking the degeneracy 

To determine flko we need oth er dy namical information. We here show that Qko can be determined if we use a new observational 
technique proposed by Loeb( 1998 ). He pointed out that the time variation of cosmic redshift might be detectable through 
two observations of ~ 10 2 quasars set a decade apart with the HIRES instrument of the Keck 10 meter telescope. The key 
point of his proposal is to use an existing spectroscopic technique, which was recently employed in planet searches. If a jupiter 
size planet exists around a certain G/F type star, the G/F type star is moving around the center of mass of the star-planet 
system with rotational velocity of the order of ~ 30m/s. In principle, the wave length of an absorption line in the spectrum of 
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the star is doppler shifted. By observing only one absorption line it is difficult to determine a period and an orbital velocity 
to obtain the mass and the orbital radius of the planet because the expected shift of the line is too small compared with the 
line width. However if one observes many absorption lines, the signal-to-noise ratio(S/N) will increase in proportion to y/N 
where N is the number of lines observed. The same technique can be applied to cosmology in principle. Each line from a 
quasar has a large line width so that it is difficult to determine the time variation of the redshift. However by observing many 
quasars and many absorption lines it might be possible to determine one parameter Qko- If such an observation is performed, 
the change of the redshift Az can be obtained as 

Az = [(1 + z)H - H(z)] At = g(z)At (23) 

where At is the time interval of two observations. A spectroscopic velocity shift is then 

Az q(z)At , . 

Av = = — , (24) 

1 + z 1 + z v ' 

which is of order ~ 1ms -1 over At ~ 10 2 years for a single quasar. Loeb showed that it is feasible to detect the cosmic signal 
(with signal-to-noise ratio of ~ 100) only over a decade with ~ 10 2 quasars. 

Since g(z) can be determined observationally, H(z) can be determined independently from Eq.(|l^) as 

H(z) = (l + z)H -g(z) (25) 

In reality it may be hard to determine g(z) for various z. However if g(z) is determined at a certain z s quite accurately, we 
can determine Sl^o as 

n fe0 = r{z s y 2 Ho 2 [l - (dr{z s )/dz) 2 ((1 + z s ) H - g {z s )) 2 ] (26) 

Even if g(z) is not determined quite accurately for any z, if the luminosity distance d^(z) — (1 + z)r(z) is obtained accurately, 
we may determine Qko by using a maximum likelihood test for many observations of g(z), since only one parameter Qko is to 
be determined. 



3.3 Application to two component model 

Finally, we note that Eq.(pl|) and Eq.(^) can be applied to the two component model such that p = pu + px, where pu is 
dust matter and px refers to the "x-component" , and that we can determine, once Qmo is known, the equation of state (or 
the effective potential) of the x-component. For example, in the case of a fluid we have 



K 2 p x {z) = 3 

K 2 px(z) = - 



(dr/dz)' 
3 



+ 



l + z) 2 - 



(dr/dzY 



H 2 Qko — H 2 Qmo 



(dr/dz) 2 



+ (1 + 2) 



d_ 

dz \ (dr/dz) 



(1 + *) - 



3r 2 



(dr/dz)' 



+ (1 + 2) 



d 



dz \ (dr/dz) 



(27) 



(28) 



Note that the cosmological constant corresponds to the px = — px case. Alternatively, in the case of a scalar field we have 



= K 2 (i+*r 



/, \ d 



1 - r 2 H 2 Q. k 
(dr/dz) 2 

1 



dz 



dz \ (dr/dz)' 



+ 



d 



2(1 + z) 2 - (1 + 2)- 



dz \ (dr/dz)'" 



k 2 V(z) = 



1 



2 [ +Z > dz \ (dr/dz) 2 



+ 



2(1 + zf- 



3r 2 



(dr/dz) 2 + 2 ( +Z) dz \ (dr/dz) 2 



(29) 



2 3 2 



(dr/dz 

Once Oa/o is measured by other observations, we can determine unambiguously the equation of state or effective potential of 
a scalar field from the luminosity distance di,(z) = (1 + z)r(z). 



4 SUMMARY 

We have explored the possibility of determining the equation of state (or the metric) directly from the observational data. 
We have shown that by combining the observation of the luminosity distance and the observation of the time variation of the 
cosmic redshift, the equation of state of the Universe (or the metric) might be directly determined unambiguously in future 
without any theoretical assumption on the matter content and the geometry of the Universe. Since there exists no convincing 
candidate for the equation of state of the x-component or the effective potential, it may be more promising to determine it 
directly from observational data. 
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